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parameteis of sensor elements in a sensor array. 
The method eomprises receiving an output signal 
of at least two sensor elements signal in reaction 
to an input signal from a signal source; estimating 
a crt^ss-conelation between the output signals 
of at least two of said sensor elements; and 
optimising a difference between the estimated 
cross-correlation and a cross-correlation model; 
and thereby estimating said parameters from 
the optimised difference. A cross-correlation 
model is used as represented by the mathematical 
equation: R=G B G"-f D 
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Title: Calibration method, device and computer program. 

rn?.nTT T>nnAT. field and BACKGROUND OF THE IlSfVENTION 

The invention relates to a method and a device for calibrating parametera of 
sensor elements in a sensor array. The invention also relates to a computer program 
product for calibrating parameters of sensor elements of a sensor array when run on a 
5 programmable device and to a sensor array calibrated with a calibration method* 

In sensor array systems, the, complex, receiver gains and sensor noise powers 
of the sensor elements in the sensor array are initially unknown , and have to be 
calibrated. (Gain) calibration enhances the quality, specifically the sensitivify, of the 
sensor system and, moreover, improves the effectiveness of array signal processing 

10 techniques for interference mitigation. 

Non-polari2sed and single polarisation calibration techniques for sensor array 
systems are generally known in the art [1,2,3,5,13], and the statistical performance 
also is a well studied topic [1,2,3,5,8,13,14]. Recently, Hamaker, Bregman and Sault 
[10, 111 12] developed, for radio astronomical purposes, a polarization formalism in 

15 which the polarization state of the received signal, and the propagation of these 
signals through the atmosphere and through the sensor array, were thoroughly and 
elegantly incorporated. This formalism is based on optics [15,16] and on extensions of 
the (approximate) solutions in radio polarimetry [17,18]. 

In tiiis formalism, the polarized signal is described by a 2x2 si2e Stokes matrix 

20 [17,18] (a Stokes matrix describes the polarization state of the signal: intensity, 
linearity, ellipticity, polarization angle, total polarizationX and the distorting and 
propagation effects by a 2x2 size Jones matrix [10,11,12] which in general is different 
for each of the dual polarization array sensors. The output of a dual polarization 
channel is described by a multiplication of Jones and Stokes matrices. The polarized 

25 array formalism is further focussed on pair-wise correlation products involving 2x2 
size Jones and Stokes matrices. However, solving systems based on this formaUsm 
require an iterative approach and convergence is not always guaranteed. Hence, in 
such systems stabiUty is a problem as well. 
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The ainele polarization and non-polarized sensor array parameter estimation 
problem is weU knovm from literature 12.13). However, these calibration methods are 
dWvantaeeoua because they require a large amount of processing. They also require 
a good initial point (gain and noise values), which ia not always available Typically 
the number of processing steps involved scales with the third power of the number of 
sensor elements. 

Recently. &st and closed form single polarization calibration techniques were 
described in [1]. In this publication. Ae calibration techniques involves the 
comparison of an estimated signal will, s signal outputted by telescopes in a telescope 
array. By optimisation of the estimated signal with a least square error 
minimalisation of the diference between «^ estimated and ouiputted signal, the 
gams of the telescopes can be derived. In the publication [IJ several variants of the 
least square error minimalisation are described. One of the minimaliaations is a 
logarithmic minimalisation. in which di&rence of the logarithms of the ct^raxiance of 

the estimated signals and the covariance of the outptttted signals are compared 

The aumber of processing steps for the logarithmic minimalisation described in 
[X] scales with the square of U« number of elements and is thus much much fester 
than conventional methods. However ttus logarithmic minimalisation has the 
^sadvantage that, for unequal gains, the method ia not efficient, which means that 
the estimation accuracy is lower than the theoretical bound. 

SUMMARY OF TFTV. TN VENTTOInT 

It is a goal of the invention to provide a better calibration me^od. more 
specific It IS a goal of the invention to provide a numericaUy stable calibration method 
by providing a method for which closed ibrm so2«to„s exist for non-polarized, single 
polariaation, and dual polarization sensors. 

Therefore, the invention provides a caUbration method according to daim 1 

Ther '""^ ''"^^^'^ -l«^ns exist 

Thereby, the parameters can be estimated in a single and stable process 

In one embodin^ent. a method according te the invention is applied to the 

cahbra^on of a dual polarised sensor array, that is an array for the reception o^d^ 

Polansed signals. For polarised signals in particular, a method according to tTe 
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invention gives new insights in the gain estimation issue, specifically, new (closed 
form) solutions to the estimation problem now become available. Furthermore, a dual 
polarization calibration method according to the invention has the advantage that 
closed form solutions exist. Such a system is numerically more stable than the known 
5 calibration systems 

In anotlier embodiment» a method according to the invention is applied to a 
single polarization or a non-polarized sensor array and a weighted logarithmic 
minimalisation is used to estimate the parameters. It can be shown, that this 
embodiment, for low S>JR is asymptotically efficient. Thus, for a large ntunber of data 

10 samples, the estimator is efiScient for low SNR. The improvement is larger for larger 
gain differences, thus the embodiment is especiaUy suited for sensor arrays with large 
tmequal gain magnitude differences. 

Furthermore, the invention provides a calibration system according to daim 
17* Also, the invention provides a calibrated sensor array according to claim 20. 

15 The invention also provides a computer program product according to claim 21. 

Such a computer program enables a programmable device to calibrate a sensor array 
in a stable manner when the program is run on the programmable device. 

Specific embodiments of the invention are set forth in the dependent claims. 
Further details, aspects and embodiments of the invention will be described with 

20 reference to the attached drawing. 

RmnBF DESCRIPTION OF THE DRAWINGS 

Fig. 1 schematically shows an example of a first embodiment of a dual 
polarization calibration system according to the invention. 
25 Fig. 2 schematically shows an example of second embodiment of a single 

polarization or non-polarized caUbration system according to the invention. 

Fig. 3 schematically shows a flowchart of an example of a calibration method 
according to the invention. 

Fig. 4 schematically shows a flowchart of an example of a calibration method 
30 according to the invention. 
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DETAILEn nKSnff,fp>n^^y 

In «»e fejlowiug description of examples of embodiments of the invention, the 
foUowing notation is used. Bold capital letters indicate matrices, bold lowercase 
letters represent vectors, and non-bold latter (either lowercase or uppercase) indicate 
scalars. Matrix and vector elements are denoted by subscripts. The subscript n 
denotes the n-th observed time sample. Superscript * denotes the compkx coniugate 
Superscript t denotes the vector or matrix transpose, which is an operation which 
switches the columns and r«ws of the matrix or vector. i.e. ((R.), = R,). Superscript H 
denotes the Hermxtian (ootflugate) transpose of a matrix ((Ea)« = ((E*% = R*«) Re{ } 
represents the real part. Im{ } the imaginary part and i is the square root of -1 E{ } 
denotes the expected value of the covariance matrix, xn denotes the output vector at 
W n. Dual polarization veete« and matrices are represented by slanted letters- for 
the single polarisation case non-slanted letters axe used. The number of single or dual 
16 polarization array sensors is denoted by p. 



6 



10 



20 



25 



Example A. Polarisation calibration 

In the following, by way of example, a gain calibration system and a method 
according to ti»e invention are described in a radio-astronomy application. However 
the invention is by no means limited to applications in radio-astronomy and may 
likewise be applied in other fields, such as for example array signal processing 
systems in which a. dominani, point like emitting source with an arbitrary 
poWation state is present, such as for example satelUte tracking phased arrays or 
otherwi^. ^ i3 based on observations of three separate or consecutive 

(polarized or non-polariaed) point-like signal sources. 

For one signal source. R represente the W or expected 2p.2p covariance 
matrix (of the 2pxX dual polarization array output vector «(t) = (x.«(t) xx^m 
xp''(t). xp^(t)). thus: ' ^ ' ^ ^' 



30 

n = E{x(t)x(t)«} ^^^^ 

.St«l. ^^"^f f ^ P-P-^- °^*he source [15.16.17.18] are represented by a 2x2 
(Stekea) matrix B. the 2px.2p noise matrix D is defined by D = mmdm, and m » 
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(diKt), diHt) , ... t dp*(t). V(t)), where dKt) and ^(t) indicate the sensor noise. D may 
be dither diagonal (in case of a perfect insulation between the two polarizations of a 
dual sensor) or block diagonal (in case of a leakage between the two polarizations of a 
dual sensor). G represents the polarization gain matrix of dimensions 2px2. 
6 Fig. I shows a dual polarisation sensor array with 2 p sensors elements by way 

of example of a system 1 according to the invention. The system 1 is aUe to perform 
the example of a method according invention illustrated in fig. 2. In fig. 1, three point- 
like sources 2i, Zss, 23, emit consecutively each a signal 3i, 82, 83 to the system 1. In 
step 101 in figure 3, each of the three consecutive signals 3i, 32, Sa is observed 

10 separately by a number of dual i>olaxi2ation sensor elements .4i,4a 4p ^ 

representing the total number of dual polarisation sensor elements) of a sensor array 
4, which may for example (but not necessarily) be an uniform linear sensor array* The 
outputs of the sensor elements 4i,43. ...»4p are connected to inputs 51 of a gain 
calibration device 6. The outputs of the sensor array are further, as indicated with the 

15 sti4>ed lines to a, not shown, beasoiformer device as is for example known in the arts 
of radars, acoustic arrays, and radio-astronomy. 

The output CTg"**^^ ^(^) contain the input signal s(t) of a single source 
multiplied with a (complex) number Ox which represents the path time delay for the 
respective sensor elements 4l The signal s(t) also multiplied with a gain factor Gi of 

20 tihe respective sensor element 4i. The output signal Xi(t) also comprises system noise 
di(t) added in each sensor element 4i to the signal s(t). 

Each dual polarization sensor 4i,42, .».,4p comprises two sensor elements, one 
for each polarization component in the signals Si-Ss. Thus> the total number of sensors 
elements in the system 1 is 2p. In steps 102 and 103> the output signals Xi^t) resp. 

25 ^{t) are presented at the outputs of the dual polarization array sensor elements 4i,42, 
..*,4p, where x^ is the output signal for one polarization component of the input signal, 
and is the output signal for the other polarization component of the input signal. 

In the described examples, it is assumed that the look direction, i.e. the 
direction of the source 2 with respect to the sensor array 4, is known. Thus, without 

30 loss of generality, the path time delay a, can be set to I, e.g. the delay is taken to be 
the same for all sensor elements. However* the invention may likewise be applied if 
the path time delay differs for some or all sensor elements. The inputs 51 are 
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communicatively coiuected to a vectoriser device S2, In step 104 the output stgnaia 
xi(t) are stacked in the output vector x(t) = (xi»(t), xi»(t) , , Xp«(t), x^it)) by the 
vectoriser device 52. 

Comiected to the vectoriser device 52 is a correlator device 53 which is able to 
S crosa-cortelate> in step 105 in fig. 3, elements in the output vector x(t). In the shown 
correlator device 53^ a covarianoe matrix is fimned by determining the covariance of 
the element ouiputa Xi(t). However other correlation methods may be used as welL 
After steps 101-105, the gain parameters, as represented by a matrix G of si^e 2px2 
(except for an arbitrary phase offiset term), are estimated in steps 106 and further, by 
10 solving cost fimctions, each cost function corresponding to a different one of the input 
signals 3i, 3%, 3s. 

The caHbration device 5 further comprises a estimator device 54 connected to 
the cross-corr^tor device 53, which is able to estimate an estimated cross-correlation 
In step 106, the estimated cross-correlation matrix IUm is estimated by the estimatar 

15 54 from the covariance matrix established by correlator device 53. However, the 
estimator device 54 may likewise estimate the cross-oorralation in a di^ront manner 
In step 107, an optimisation device 55 cdmmiuucatively connected to the 
estimator device 54, compares the estimated covariance matrix JU^f with a model of 
the covariance matrix R stored in a memory 56. 

20 In the example, in step 107 a difference between the cross-correlation model 

and the estimated cross-correlation matrix Ret^ is optimised, by minimising • the 
difference, which difference is also referred to as the estimation error. In step 108 the 
gain parameters are estimated or derived from the result of step 107. In the example 
of fig. 3, in step 106 the gain parameters are derived from the optimised difference 

2S between tlie estimated matrix and the model by a gain estimator device 57. The 
estimator device 57 is connected to the optimisation device 55, The estimator device 
57 is also connected to an calibrator output 58, via which the estimated parameters 
can be transmitted further to be used in further processing such as the processing of 
signals -received at the sensor elements after the calibration. 

30 the example of a device of fig. 1 and the example of a method of fig. 3, as a 

crofis^correlation model, a model of the covariance matrix -R is used. The covariance 
matrix Ri for one source can be expressed as: 
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Rx=GBiG^ + D (A2) 
and for Q sources the covariance malacix 22 can be expressed as: 

ie= {Sm^»..QGB»GH}+J). (A3) 

Using the model of equataons A2 and/or A3, given an estimate covariance 
5 matrix Beat, for one or more consecutive sources, the gain matrix O and noise matrix 
D can be estimated. In this example, in step 106 an estimate of the covariance matrix 
is obtained by the estimator device 54 by forming the product: 

JJ»rf= 2n X (t n) X(tn) « (A4) 

Thus, in this example A, the estimated covariance matrix i!nt is a time 
10 averaged covariance matrix of the output signals. However, it is likewise possible to 
estimate tiie covariance matrix in a different manner. 

It can be shown [12J that for a system of three di£Eerent]y polarized sources an 
unambiguous or unique gain solution can be obtained (apart firom an arbitrary phase 
o£&et valid for the entire system). The polarized gain can be found by solving in steps 
15 107 and 108 the equation: 

{GM = argrain oj, {(^Bc^i - {GBiG^D >0p)''+ 
(||j2«.2-{£3lB2GJM-2>}ir)» + (!|J8e«.s-{GBsGH+i>>||j.)s) (A6) 

20 Here the subscript P means the Frobenius norm, which ia tiie square root of 

the sum of the squares of all matrix elements. Thus, the added least square 
differences between the estimated covariance matrices for output signals resulting 
from input signals from three sources and the corresponding covariance matrix 
models are optimised. However, it is likewise possible to optimise the difference in a 

25 different manner. The estimated gain and noise parameters of the sensor elements 
are then estimated to be the gain and noise parameters for which equation (A5) is 
solved. 

As an example of the solution of eq. (AS), a rank two factor analysis approach 
may be ^xaed. as is briefly described below. However, other more or less standard 
30 solutions for the equation (A5) exist and the invention is not limited to the example 
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described below. Solutions for equations similar to eq. <A5> are generally known in the 
art. The example method starts with solving: 



Ri^GBiG^^D (AS) 
5 for the estimated covarianoe matrix iSest^i of the first signal. 

This can be carried out by applying a rank two fkctor analysis [6] to 
submatrices M of Itea^i (both M and liest^i are rank two), in which the main (block) 
diagonal of R^i is not included as it contains the system noise as weU. A resulting 
(rank 2) matrix A can be formed, based on the obtained Jfcf s» with the properly that, 
10 for the off (block) diagonal elements^ AA^ is equal to R, D can be obtained by 
calculating Rm,i - AA^. 

The next step is to calculate the matrix root of B so that equation [A6] can be 
expressed as: 

R « (GB^ (GB^^ H + Zl (AT) 

X5 

With the factor analysis approach, described above, a GB^ foctor can be found. 
The problem here is that the solution is not imique as any unitary matrix can be 
inserted in equation A7 (between the two GB^ factors) without affecting B. It can be 
shown that equations like equations A6 and A7 can be fully solved (up to a single 
20 arbitrary phase offset) by using three subsequent, independent, observations (IZast.!, 
iZftct.s, and Re^^) of emitting point soiurces with different polarisation states using a 
geuerali2sed eigenvalue analysis and the pseudo inverse. 



Example B: Gain calibration by weighted logarithmic estimation 
25 In fig. 2, an example of a system 1 according to the invention is shown. In fig, 
2, a point-like source 2 emits a signal 3 to the system 1* The system 1 comprises 
sensor elements A\M> »*-»4p (p representing the total number of sensor elements) of a 
sensor array 4, for example (but not necessarily) a uniform linear array. The sensor 
elements receive an input signal 8(t) which stems from the signal 3. At the outputs of 
30 each of the array sensor elements 4i,42, .•.,4p an output signal xi(t) is presented. The 
outputs of the sensor elements 4i,42 4p are connected to a gain calibration device 5. 
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The output signals aa^t) contain the input signal si(t) multiplied with a 
(complex) Dxunber Sr which represents the path time delay for the respective sensor 
element 4i The signal Si(t) is also multiplied with a gain factor gi of the respective 
sensor element 4u The output signal xi(t) also comprises syBtem noise di(t) added in 
5 each sensor element 4i to the signal 8i(t). 

The gain calibration device 5 is able to perform the example of a method 
according to the invention of fig. 4, The calibration device 6 has a number of inputs 51 
connected to the outputs of the sensor elements 4i,43, -..,4j?- In step 201, the output 
signals 7ci(t) are received ficom the sensor elements. The inputs are commuxiicatively to 

10 a vectoriser device 52 in which in step 202 the output signals Xi(t) are stacked in a 
vector, denoted by bold lowercase letters: x(t) = [xi(t), ... , Xp(t)]^ 

In the fbUowii^, it ia assumed that the look direction, Le. the direction of the 
source 2 with respect to the sensor array 4» is known. Thus, without loss of generality 
the path time delay dt can be set to 1, e.g. the delay is taken to be the same for all 

15 sensor elements. However, the invention may likewise be applied if the path time 
delay differs fbr some or all sensor elements. The array ou^ut vector s:(t), under the 
assumption mentioned above, can be described as x(t> = e ^(t) + d(tX where d(t) is a 
vector containing the array noise signals di,d2»...,dp and where g denotes the 
(complex) array gain factors gi,g2,...,g|> of the sensor elements 4i>42, .•.,4p- 

20 Cozmected to the vectoriser device 52 is a correlator device S3 which is able to 

cross-correlate, as in step 203 in fig. 4^ the elements in the output vector ac(t). In the 
shown correlator device 53 a covariance matrix is formed by determining the 
covariance of the element outputs xi(t). However other cross-correlation methods may 
be used as well. The true covariance matrix R is the expected value of the covariance 

25 matrix, thus: 

B^E{x(t)x(t)H} (Bl) 

The calibration device 5 further comprises a estimator device 54, which is able 
30 to estimate an estimated cross-correlation, as in step 204 in fig. 4. In the shown 
example, the estimated covariance matrix Rest is estimated by taking the, weighted, 
time-average of the covariance matrix of the output signals Xi(t) from the correlator 
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device 6S. However, the covarianoe matrix can likewise be estimated in a di£ki«nt 
manner. In the shown example, the estimated ooyariance matrix R«<t can 
mathematically be desoibed as; 

R«t= (1/N) En x(t«) x(t») H (B2) 

In which equation (B2), N represents the total niunber of the ouipiit vector 
x(t). Since the noise is imcorralated to the source signal, the eovariauce matrix R can 
be modelled, see reference [1], as: 

B»SgB + D (B8) 



Here, the noise contribution is represented by a diagonal matrix D with the 
noise variance on the main diagonal ( D ^ £»{ d ). Thiis. ficom the estimated 
15 oovariance matrix lUnt. an estimated gain gMt and noise Dcm can be deduced, because 
6«se=» e»st gest K-f Dest. In tilis model |;i] it is assumed that the source power is unity. 

The estimated gain and estimated noise are then derived by ininiTn^fi ine in step 
20S and 206 by an optimisation device 55 commimicatively connected to the estimator 
device 54 and a memory device 56 in which the model is stored, a difference between 
20 the estimated oovariance matrix R«m and the model of the covariance matrix B, which 
difference is also referred to as the estimation error. A generally used approach 
thereto is to minimize the least squares (LS) cost function: 

igest,I>^} = argmin«i> ( | Re,t - (g g» + D) j f) * (B4) 

25 

where subscript P denotes the Frobenius matrix norm. In [1] the minimization is done 
by solving 

feast) = argminsJi; ( J J vec(ln(Re*t) - ln(g gH) + 2xki) || p)2 
where J is a selection matrix which pute zeros on the main diagonal elements, 
30 thereby removing the system noise contributions, where k is a phase unwrapping 
vector containing integer values, and where vecO is a matrix operator stacking the 
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matrix elements in a vector. Once the gains are found, the noise matrix is found 
by 

D«<t = Bert-g6rtgert» CB6) 

5 In the example of figs. 2 and 4, a weighted logarithmic least square cost 

ftmction (WLOGLS) ia used. The WLOGLS is equal to: 

{ges,) = argminB*(|| WJvec(hi(R«0-hifegH) + 2wki) ||r)« CB6) 

10 W represents a weighting matrix equal to a)-»« n®a>-'« V). in which r is 

defined as diag(lg»l....|gpl). ie- T is a diagonal matrix with values of the diagonal 
elements set to the absolute value of the, estimated, gain factors gi. g2..... g*. Hei« ® 
denotes the Kronecker matrix product. For a low Signal to Noise Ratio (SNR), the 
SNR is substantially equal to g g«/ trace(D), it can be proved (not shown here) that 

15 the gain estimation is asymptotically efficient, which means that the estimation 
accuracy, for a large number of observed samples, meets the theoretical bound (i.e. it 
means that it theoretically can not be estimated more accurately) 

To obtain a fest and reliable estimate, the least square optimisation can be 
performed in a number of times. In the shown example, the optimisation is performed 

20 two times in steps 206 and 206. In step 205, an estimate is obtained of the gain g and 
noise D by observing the point source and using the identity matrix I as a weight 
matrix. In step 206, the gains are estimated using the weight matrix constructed as 
described above using the estimated gain ge«t and noise D«sc estimated in step 205. 

The calibration method of this example B works specially weU for low SNR 

25 and als^) for arrays in which the antenna gains have a wide gain magnitude variation. 
The latter ii also useful ft* gain estimation in the case part of the antennas (due to for 
example malfunctioning) have very low SNR. 

A method or device according to the invention may be applied in any array 
signal prpcessing system, such as systems in which a dominant point-like emitting 

30 source is present. For example, a method or a device according to the invention can be 
used in cellular telephone base stations, phased array antennas or otherwise. Also, 
the invention may be applied for calibration of a directional hearing device. 
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comprising an array of microphDnes which can for example be phase-tuned to receive 
acoustical signals ftom a certain direction and reject or filter out signals from other 
directions. Furthermore, the signals to be sensed by the sensor array may be of any 
suitable type, such as &r example radio signals^ acoustical signals^ optical signals or 
5 otherwise. Likewise, the signal source may be of any suitable type, such as a satellite 
in orbit around a celestial body or a pulsar. In general , pulsars are extremely precise 
pulsatixig celestial bodies and hence with a sensor array directed to a pulsar a precise 
tlme-measiutog device can be obtained. 



10 comprising data representing a computer program product, comprising program code 
for performing steps of a method according to the invention when run on a 
programmable device. Such a data carrier can for example be a read only memory 
compact disk or a signal transfer medium^ such as a telephone cable or a wireless 
connection. The programmable device may be of any suitable ^e. For example, it 

15 may be a computer communicatively connected to a senor array. However, the 
computer may likewise be not connected to a sensor array, but receive data 
representing s^nals from the array, e.g. via a floppy disk or a compact disk. 

It shoiild be noted that the above-mentioned embodiments illustrate rather 
than limit the invention, and that those skilled in the art will be able to design may 

20 alternatives without departing firom the scope of the appended claims. In the claims, 
any reference signs placed between parentheses shall not be construed as limiting the 
claim. The word 'comprising* does not exclude the presence of other elements or steps 
than those listed in a claim. The mere fact that certain measm^es are recited in 
mutually different claims does not indicate that a combination of these measures 

25 cannot be used to advantage. Furthermore, if a document is referenced to in this 
application, this does not indicate that the document relates to the same field of 
technology as the present invention. 



Furthermore, the invention can likewise be applied as a data carrier 
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CLAIMS 

1. A method for calibrating parametera of sensor eiemente in a sensor array, 
comprising: 

5 receiving an output signal of at least two sensor elements signal in reaction to an 
input signal from a signal source; 

estimating a cross-correlation between the output signals of at least two of said sensor 
elements; 

optimising a difference between the estimated cross-correlation and a cross- 
10 correlation model; and thereby estimating said parameters from the optimi^d 
difference; 

wherein a cross-correlation model is used as represented by the mathematical 
equation: 

15 in which equation: 

R represents a cross-correlation matrix, 

6 represent a gain matrix comprising gain parameters, 
represents the Hermitian conjugate of the gain matrix, 

D represents a ((block)- diagonal noise matrix comprising noise parameters and 
20 B represents a matrix comprising information about the signal source. 

2. A method as claimed in claim 1, wherein said difference is a least square 
difference. 

25 3, A method as claimed in claim 1 or 2, wherein the cross-correlation is obtained 
by determining a time-averaged covariance matrix from the output signals. 

4. A method as claimed in any one of the preceding claims, wherein the sensor 
array is a single polarisation or non-polarized sensor array. 

30 

5. A method as claimed in claim, wherein the sensor elements are dual 
polarization sensor elements for receiving a dual polarised signal 
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6. A method as claimed in any one of the preceding claims, wherein said method 
is performed for output aignals of the sensor elements generated in reaction to input 
signals firom at least three signal sonroes with different polarizations. 



5 7. A method as claimed in claim 4, wherein said optimising comprises: 
minimising a difference between a weighted logarithm of the estimated cross- 
oorrelation and a weighted logarithm of the cross-correlation and 
estimating the gain of at least one of the sensor elements firom said difference. 

10 8. A method as claimed in claim 7, wherein the logarithm is weighted by a 
weighting matrix with matrix values relating to said gain parameters! 

9. A method as claimed in claims 7 or 8, wherein said optimising and said 
estimatix^ gain parameters are performed at least a first time and a second time, 

15 wherein in the first time an uniform weight is used for all output signals and in the 
second time the weight is used in dependence on the gain estimated in the first time 
for the respective output signals. 

10. A method as claimed in any one of claims 7-11, wherein said optimising 
20 comprises an operation as represented by the mathematical equation: 

{g«t} = argmin«jr (II WJ vec(ln(B«st) -ln(ggH) + 27rki) 
, in which equation: 

geAt represents the parameter to be estimated; 
25 g represents a variable 

g« represents the Hermitian conjugate of the variable 

J represent a selection matrix which puts zeros on the main diagonal; 

k represents a phase unwrapping vector containing integer values. 

W repre^^ents a weighting matrix; and 
30 'Rest represents the estimated cross-correlation. 
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11. A method as claimed in any one of the preceding daims, wherein the signal 
source is a satellite in orbit around a oeleatial body. 

12. A method as claimed in any one of the preceding claims, wherein the signal 
source is a pulsar. 

13. A method as claimed in any one of the preceding claims, wherein the output 
signals have a low signal to noise ratio. 

14. A method as claimed in any one of the preceding claims, wherein the sensor 
elements are antennas in a phased array antexma. 

15. A method as claimed in any one of the preceding claims, wherein the sensor 
elements are electro-magnetic sensors elements. 

16. A method as claimed in any one of the preceding daimsl, wherein the sensor 
elements are acoustical sensor elements. 



17. A calibration system fcr cahbrating parameters of sensor elements in a sensor 
20 array , comprising 

at least two inputs, each conneetable to an output of an sensor element in a sensor 

array; 

a correlation estimator device for estimating a correlation between the output signals 
of at least two of said sensor elements 
25 an optimiser device for optimising a diflference between the estimated cross- 
correlation and a cross-correlation model and thereby estimating said parameters 
from the optimised difference; 

a memory device containing the cross-correlation model, which model is represented 
by the mathematical equation: 
30 ifc=GJSGM + /) 

in which equation: 

R represents a cross-correlation matrix. 
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G repiresent a gain matrix comprising gain parameters, 

represents the Hermitian conjugate of the gain matrix, 
2) represents a noise matrix comprising noise parameters and 
B represents a matnx comprising information about the signal source 
and , 



18. A calibration system as claimed in claim 17, wherein the sensor array is a dual 
polarised sensor array 

10 19. A calibration system as claimed in claim 17, wherein the sensor array is a 
single polarization or non-polarized sensor array, 

20- An array signal processing system calibrated with a method as claimed in any 
oneof claims 1-17 

IS 

21. A computer program product, comprising program code for per&rming steps of 
a method as claimed in any one of claims 1-17 wheii run on a programmable device. 

22. A data carrier comprising data representing a computer program product as 
20 claimed in claim 21. 
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